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We study pulse propagation in one-dimensional chains of spherical granules decorated with small 
randomly-sized granules placed between bigger monodisperse ones. Such "designer chains" are of 
interest in efforts to control the behavior of the pulse so as to optimize its propagation or attentua- 
tion, depending on the desired application. We show that a recently proposed effective description 
of simple decorated chains can be extended to predict pulse properties in chains decorated with 
small granules of randomly chosen radii. Furthermore, we also show that the binary collision ap- 
proximation can again be used to provide analytic results for this system. 
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I. INTRODUCTION 



The study of pulse propagation in granular media has 
become a field of intense research interest. This is partly 
due to the fundamental importance of understanding the 
associated nonlinear dynamics and partly due to its di- 
rect application to our day-to-day lives. Any system in 
which interactions among its discrete constituents are 
solely described by their macroscopic geometrical shapes 
and elasticity properties rather than by their microscopic 
(atomic or molecular) nature can be classified as granu- 
lar system. Among granular systems, a class that draws 
special attention is that of so-called dry granular sys- 
tems. A defining property of these systems is that the in- 
tergranule interaction potentials are always positive i.e., 
they are purely mutually repulsive. Furthermore, they 
are only nonzero as long as the two bodies are in physi- 
cal contact. This peculiarity, together with the spatially 
discrete nature of granular systems, gives rise to a host 
of interesting phenomena. Depending on parameter val- 
ues, granular systems can express liquid-like, solid-like 
(glasses) , or gas-like properties [l], Q . 

It is well known that an initial kinetic energy impulse 
imparted to an edge granule of a chain in the absence 
of precompression can result in solitary waves propagat- 
ing through the medium [j, [I|. In the recent past, pulse 
propagation in one-dimensional (ID) chains of granules 
has been studied extensively both theoretically and ex- 
perimentally [5J-[23j] - These studies arc in part inspired 
by a number of practical applications, e.g., in the design 
of shock absorbers J2J, [25( , sound scramblers [26|, [27| and 
actuating devices |28j . 

Structural variations along a granular chain can sig- 
nificantly influence pulse propagation through the sys- 



tem. These structural variations (or polydispersity) are 
frequently introduced in a regular fashion such as in ta- 
pered chains, in which the size and/or mass of succes- 
sive granules systematically decreases or increases. A 
detailed s tudy of these effects has been carried out nu- 
merically jUT - tla . |29| . In our own recent work [3(J we 
introduced a binary collision model to derive fairly sim- 
ple analytic results to describe pulse propagation in var- 
ious ID tapered chains. We showed that most of the 
essential properties that characterize pulse propagation 
in tapered chains can be extremely accurately described 
in these systems using the binary collision model. 

Recently we turned our attention to more complex 
chain configurations, specifically, to decorated chains, 
that is, to chains in which large and small granules al- 
ternate in some regular fashion. These chains can not 
be studied directly using any binary collision model, and 
the reason is quite obvious: as a pulse propagates, the 
small granules rattle back and forth between their larger 
neighbors and thus at least three granules rather than 
two are involved in elementary collision events. How- 
ever, we introduced an effective description |3l| whereby 
we represented decorated chains by associated undeco- 
rated ones. We tested this methodology on a variety of 
simple and tapered decorated chains and showed that the 
effective undecorated chains reproduced the behavior of 
pulse propagation in the original decorated ones in all 
cases provided the small decorating granules were suf- 
ficiently small (see below). This effective representation 
can then be treated analytically using the binary collision 
approximation. The resulting analytic expressions allow 
us to explore regimes in which numerical algorithms may 
be unstable, and even regimes where the pulse amplitude 
becomes so weak as to be close to the numerical noise. 

All the studies described above have dealt with chains 



with regular configurations, be they monodispcrsc or ta- 
pered, simple or decorated. Here we introduce a new 
element, namely randomness. It is interesting to study 
the effects of randomness because (I) most granular sys- 
tems in nature are not regular, and (2) randomness might 
be used as a control element in manipulating pulse prop- 
agation properties in designed systems. It therefore be- 
comes important to understand effects of randomness on 
pulse propagation. Except for recent work by Fraternali 
ct al. [32j , we are not aware of any other work on pulse 
propagation in random granular chains. In this work we 
develop an understanding of some effects of randomness 
on pulse propagation. For this purpose, we extend our 
previous work on decorated chains [31( to incorporate ef- 
fects of randomness. 

Our randomly decorated chain is constructed from a 
monodispcrsc array of large granules separated pairwisc 
by smaller granules whose size is randomly picked from 
a pre-assigned distribution. For our effective descrip- 
tion [31j to be valid, the radius of the smaller granules 
must be no larger than ~ 40% of the larger granules. This 
places a restriction on the upper cutoff of the size distri- 
bution of smaller granules. We show that a randomly 
decorated chain can then be mapped onto an effective 
undecorated chain with random masses and random in- 
teractions. We find that for all properties studied herein, 
the behavior of our effective chain (obtained from the nu- 
merical solution of the equations of motion for this chain) 
is in remarkable agreement with the behavior of the orig- 
inal decorated chain (obtained from the numerical solu- 
tion of its equations of motion). In particular, we focus 
on three properties of the pulse. First, we compute the 
pulse amplitude and its variation along the chain. Sec- 
ond, we determine the average speed of the pulse along 
the chain. Third, we calculate the distribution of the 
times that the pulse takes to reach the end of the chain. 

We then go on to test our binary collision approxima- 
tion applied to the effective chain. We find that whereas 
the approximation does not reproduce the pulse ampli- 
tude well (for reasons that we understand and for which 
we suggest a possible remedy), the other two properties 
are extremely well predicted analytically by the model. 
This provides a powerful tool to avoid costly numerical 
simulations. 

In Sec. |TI] we introduce the granular chain model in 
terms of rescaled (dimensionlcss) variables. In Sec. IIIII 
we introduce decorated chains and present our effective 
description in terms of undecorated chains with renor- 
malized interactions and masses. The analysis here is a 
generalization of the work of Ref. [3l| to accommodate 
the random variation (within limits) of the radii of the 
smaller granules. In this section we also exhibit the an- 
alytic results obtained by applying the binary collision 
approximation to the effective chain. Comparisons with 
numerical results are presented in Sec. IIVI In Sec. IV1 we 
provide a summarizing closure. 



II. THE MODEL 

We consider chains of granules all made of the same 
material of density p. When neighboring granules col- 
lide, they repel each other according to the power law 
potential 



V = -r' k \y k -y k+1 \ n . 
n 



(1) 



Here y k is the displacement of granule k from its posi- 
tion at the beginning of the collision, and a is a con- 
stant determined by Young's modulus and Poisson's ra- 
tio [33|, |34[. The exponent n is 5/2 for spheres (Hertz 
potential), which we use throughout this paper in our 
explicit calculations [34J|. We have defined 
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where R' k is the principal radius of curvature of the sur- 
face of granule k at the point of contact. When the gran- 
ules do not overlap, the potential is zero. The equation 
of motion for the A:th granule is 



M, 



d 2 y k 
dr 2 



a? i-i(yfc-i -Vk) n 1 d(y k - 



Vk) 



ar' k (y k - Vk+i)" l 6(y k -y k+ i), (3) 



where M k = (4/3)7r / o(i?J c ) 3 . The Heaviside function 9(y) 
ensures that the elastic interaction between grains is only 
nonzero if they are in contact. Initially the granules are 
placed along a line so that they just touch their neighbors 
in their equilibrium positions (no precompression), and 
all but the leftmost particle are at rest. The initial ve- 
locity of the leftmost particle (k — 1) is V\ (the impulse). 
We define the dimensionless quantity 
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and the rescaled quantities x k , t, m kl and R k via the 
relations 
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M k = Mim k . 



y k = R'^^Xk, 

R k = R 1 Rk, 

Equation ([3]) can then be rewritten as 

rk-i(xk-i - x k ) n ~ 1 9{x k -i - x k ) 
-r k (x k - x k+1 ) n ~ 1 d(x k - Xk+i), 



m k x k 



(5) 
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where a dot denotes a derivative with respect to i, and 
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The rescaled initial velocity is unity, i.e., «i(t = 0) = 1. 
The velocity of the k-th granule in unsealed variables is 
simply V\ times its velocity in the scaled variables. 



III. THEORETICAL RESULTS: EFFECTIVE 

CHAIN AND BINARY COLLISION 

APPROXIMATION 

A. Effective chain 

We consider a decorated chain with a small granule in- 
serted between each pair of large granules. The end gran- 
ules are large, and the large granules are monodispcrse. 
The sizes of the small granules arc random. In order 
to obtain an effective description for the pulse dynamics 
in the decorated chain, wc follow the scheme presented 
in Ref. [31(. In that work we showed that the effective 
description could be constructed using the outcome of a 
detailed analysis of a chain of five granules. This five- 
granule chain provides the elements of all the granules 
in the actual long chain: large granules in the interior of 
the chain, small granules in the interior of the chain, and 
large granules at the ends of the chain. 

Consider, then, a chain of five granules labeled from 
k — 2 to k + 2 in unit steps, granules k—1 and k + 1 being 
the small granules. The radius of the large granules is 
Ri = R = 1 in rescaled variables, and those of the small 
granules is Rk = R' k /R[. The dynamics of this chain of 
granules is governed by the set of equations, 

-r fc _ 2 (x fe _ 2 - x k -i) n ~ 1 9(x k - 2 - Xk-i) 
r»_i(a:;_i - Xi) n ~ 1 6(x i -i - Xi) 
n(xi - x i+ i) n ~ 1 8(x l - x i+ i), i — k, k ± 1, 
rk+i(xk+i - x k +2) n ~ 1 0{:xk+i - Xk+i)- (8) 

Since all the large granules are of the same size, it follows 
that m,k-i = TUk+i = rnk+2 = m (= 1 in the rescaled 
units), and we have rk-\ = rk-2 and r k = rk+i, which 
just says that for a given small granule the left and right 
sides are geometrically identical. Following Ref. |3l|, this 
chain can be mapped onto an effective chain of three 
large granules which we call "left" (I), "right" (r), and 
"between" (6), with renormalized masses fii, /i r and fib 
given by 
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Hi = m+-mfc_i, fj, r = m+-m k+ i 
fib = m+ -(wfc-i + mk+i). 



(9) 



The effective interaction between two neighboring large 
granules in this three-grain configuration is 
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Now we return to the actual long chain, which is made 
of N large granules and N — 1 small ones. Our effective 
chain then has N effective granules where each of the 



two edge granules corresponds to either the "left" or the 
"right" effective granule described above, and the rest of 
the granules in the chain are described according to the 
"between" granule prescription. Wc emphasize that since 
the size Rk of the smaller granules is a random variable, 
the renormalized masses, Eq. (|9]), and the effective inter- 
actions, Eq. (|10l) . in the effective chain are also random 
variables. This completes the mapping onto an effective 
chain. This mapping in turn allows us to implement the 
binary collision approximation, which we do next. 



B. Binary collision approximation 

Using the binary collision approximation, in our pre- 
vious studies we successfully calculated the time that it 
takes the pulse to reach the end of the chain. For this 
purpose wc first calculate the time spent by the pulse 
at each granule. If we assume that in the effective chain 
the pulse propagates through a series of successive binary 
collisions, the time taken by the pulse to go from the kth 
granule to the (k + l)st granule in the chain is given by 
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and Mk,k+i = fJ-kfJ-k+i/ifi-k + Mfc+i) is the reduced mass 
for the pair of granules k and (k + 1) |3l|. The velocity 
amplitude Vk in Eq. (|12[) is 



Vk 



fc-i 
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The time taken by the pulse to reach the end of the chain 
is therefore 

N-l 

In = 2_^ Tk,k+i 
fc=i 
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Suppose now that the small granules are chosen ran- 
domly from a given distribution, i.e., their sizes and con- 
sequently their masses are random variables. This means 
that each /c-dcpcndcnt term on the right hand side of 
Eq. (|12j) is a random variable. To help us manage these 
contributions, we implement two further approximations 
for Eq. (fT5"j) . First, since the radius of the small granules 
is restricted to be no larger than 40% of that of the large 
granules, Rk < 0.4, we neglect the correction of order R^ 
to the masses of the (large) granules in the effective chain. 
All granules in the effective chain then have equal (fixed) 



masses equal to that of the large masses in the original 
chain [see Eq.©]. Thus Vk = 1 and M.k,h+i = m/2. 
Equation (fT5|) then reduces to 



t N = A(n) (-J £ 
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(16) 



Secondly, we assume that Rk is sufficiently small to ne- 




FIG. 1: (Color online) Averaged velocity profile of the larger 
granules with N — 11. Results from the effective scheme are 
shown as black continuous curves. 

gleet it in the numerator of the above expression, which 
is consistent with the requirement for the validity of the 
effective description in the first place. We thus approxi- 
mate Eq. ([To]) as 



t N « A(n) 
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(this assumption might be questionable when Rk is as 
large as 0.4, but the results shown later support the ap- 
plicability of this approximation). 

Although all the formulas given above are valid for any 
distribution function subject to the stated restrictions, 
for concreteness we implement a uniform distribution of 
small radii over the interval [a, 0]. It can easily be shown 
that each term Xk = (Rk)^ 1 ^ 2n on the right hand side of 
Eq. (|17JI is then distributed over the range [f3~~ ,a~"^] 
according to 



P(X k ) = 



2n 



-X, 



-(2n+l) 



(18) 



Since this distribution has finite first and second mo- 
ments, as TV — > oo it follows from the central limit theo- 
rem that the sum over the independent random variables 
Xk tends to a Gaussian distribution. The mean and the 
variance of the Gaussian are obtained from the mean 
(k) and the variance (a 2 ) of the underlying distribution. 
Thus for large N, t^ is distributed around the mean Nk 
with a variance given by Na 2 . It is straightforward to 
see that the average time (tjv) then varies linearly with 
N, (iiv) ~ S(N — 1), and the slope of the line is given by 
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(19) 



We stress that Eq. (fl7|) and consequently (|19p are valid 
when the effect of randomness is negligible on the masses 
of the large granules and is only important in modifying 
the effective interactions between them. Note that within 
the binary collision approximation these expressions can 
be used not only to calculate the time that the pulse 
takes to arrive at the end of the chain but to arrive at 
any specified granule of the chain. 

In the next section we shall compare these results with 
the numerical solution of the exact equations and also 
with the numerical solution of the equations for the ef- 
fective chain. 



IV. COMPARISON WITH NUMERICAL 
RESULTS 

To restate our scenario, we consider a chain of N iden- 
tical large granules, each pair of which is separated by 
a small granule. In all our chains the first and the last 
granules are of the large variety. The radii of the (N — 1) 
small granules (or masses) are randomly selected from a 
uniform distribution in the range [a,/3]. For all of our 
numerical results we choose a = 0.01 and P = 0.4. The 
value of P = 0.4 (and no larger) is dictated by the va- 
lidity of the effective description. Note that for = 0.4, 
the mass of granules in the effective chain (Eq. [9]) can 
randomly increase by up to ~ 6% of the mass of the 
large granules in the original chain. The density of all 
granules is the same. The equation of motion for each 
granule is given by Eq. ©. We solve the set of 2N — 1 
coupled equations numerically for a large number of re- 
alizations (~ 30000) of the radii of the small granules. 
In this section we compare these solutions with the cor- 
responding results for the effective chain, and also with 
those obtained from the binary collision approximation 
to the effective chain. 

In Fig. Q] we show the velocity amplitude profile of the 
bigger granules averaged over all realizations. We observe 
that a well-behaved (average) pulse of time-varying am- 
plitude and width propagates through the random chain. 
In the same figure we show the results obtained from the 
effective description, where the granular masses and their 
interactions are given by Eqs. © and (fTU)) , respectively. 
The two results are almost identical, thus showing that 
a description in terms of the effective chain is valid even 
in the presence of randomness. 

In Fig. [5] we show the change in the pulse ampli- 
tude as it passes from one large granule to the next 
along the chain. The results are indistinguishable be- 
tween the original decorated chain and the effective chain. 
The amplitude follows a stretched exponential decay 
Vk = aexp(— bk c ) with fitted parameter values a = 0.756, 
b = 0.067, and c = 0.714. The inset in the figure shows 
the stretched exponential decay in the pulse amplitude as 
a function of time, v(t) = a'cxp(— b't c ) with a' = 0.744, 
b' = 0.035, and c! = 0.714. Note that c = c', indicating 
that time t and granule number k are linearly related, 
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FIG. 2: (Color online) Decay in the average pulse velocity 
amplitude as a function of the (larger) granular number k in 
the effective chain with N = 26. The continuous curve repre- 
sents the stretched exponential decay vt = aexp(— bk c ) with 
a — 0.756, b = 0.067 and c = 0.714. Inset: stretched expo- 
nential decay in the pulse velocity amplitude as a function of 
time t. 



i.e., t (x k. 

The variation of t with k in the original chain or in 
the effective chain (they are again indistinguishable) is 
shown in Fig. [3] for N = 26. The linear relation between 
k and t was obtained in our earlier work onpulse propa- 
gation in monodisperse granular chains 0, Q . The slope 
of this linear variation in time as obtained from the fit 
to the numerical data is ~ 2.402. This is within 2% of 
the value S = 2.348 obtained from the binary colision 
approximation, Eq. (fll?]) . Thus the effective chain as well 
as the binary collision approximation to it yield excellent 
results in agreement with those of the original chain for 
the variation of the pulse propagation time as a function 
of granule number. 




FIG. 3: (Color online) Average pulse propagation time as a 
function of k for N = 26 for the effective chain. The straight 
line is the fit t = 2.402/fc - 1.181. 

The effective chain captures the behavior of the origi- 
nal chain extremely well for all the properties considered 
above. Indeed, the behavior of the two is essentially in- 
distinguishable. We now go on to assess the validity of 
the binary collision approximation for the effective chain. 
Above we showed that the average pulse propagation 
time as a function of granule number is captured very 
accurately by the approximation. Next we focus on the 
distribution of times that the pulse takes to reach the end 
of the chain. Because of the randomness in the chain, 



this time is distributed around a mean. We compute 
the distribution of this time by solving the exact dynam- 
ical equations of motion for the decorated chain over a 
large number of realizations (~ 30000) of the distribution 
of smaller granules. The theoretical prediction for this 
time obtained from the binary collision model is given 
by Eq. (fT51) . In Fig. @] we show comparisons between the 
theory (filled circles) and the numerical results (empty 
circles) for various lengths of the chain. In showing the 
comparison in Fig. 21 we have adjusted the peak position 
of the distribution obtained from the theory, which gives 
slightly higher values of the peak position. This shift 
arises from the small error in the prediction of the pulse 
velocity when using the binary collision approximation. 
In calculating the distribution of arrival times, this small 
error is accumulated over all the terms in the sum, that 
is, it is in effect multiplied by N, the length of the chain. 
However, we note that the amplitude and the shape of 
the distribution are very well reproduced by the theory. 
For a short chain (N = 6, 11) the distribution of ar- 
rival times of the pulse at the edge of the chain is quite 
asymmetric. This asymmetry decreases as the length of 
the chain is increased, and for N = 26 the distribution is 
approaching a Gaussian. This is the result of the central 
limit theorem. Remember that in Eq. (|15[) we are adding 
random terms. However these terms are not indepen- 
dent (each term contains information about the random 
size/mass of all the previous granules through Vk), and 
therefore the sum is more complicated than one involv- 
ing independent random variables. In any case, since for 
large N the sum assumes a Gaussian-like form, the corre- 
lation between different terms in Eq. (|15|) is presumably 
small or highly localized (falling off quickly with increas- 
ing k). 
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FIG. 4: The distribution of arrival times of the pulse to gran- 
ule k = 2N — 3. This is the large granule before the last 
large granule at the end of a chain of 2N — 1 granules. Here 
N = 6, 11, 16, 26 starting from the left upper panel and chang- 
ing in the clockwise direction. Filled circles: theory; empty 
circles: numerics. 

In order to quantify the difference between the 



computed distribution functions in Fig. [4] and the 
standard normal distribution, we have performed the 
Kolomogorov-Smirnov [35| statistical test on the data ob- 
tained numerically and used to generate the distribution 
shown in the figure. For this purpose, we first collect the 
data {X n = xi, £2, • • • x n } computed from the numerical 
solution in increasing order x\ < X2 < ■ ■ ■ x n . Then the 
mean ji n and the standard deviation a n are computed as 
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m(«) = -Y] 

a(n) 



\ n — 1 z -—' 



(20) 



A new data set z n = (X n — /i(n))/a(n) is then gener- 
ated. The Kolmogorov-Smirnov test involves computing 
the statistics of the absolute difference (non-directional 
hypothesis) between the cumulative frequency distribu- 
tion F z (n) of z n and that of the standard normal distri- 
bution F , i.e., d = \F z (n) — F \ [35|. Acceptance of the 
null hypothesis that the distribution is Gaussian within a 
given level of confidence requires that d be appropriately 
small. Consistent with the progression seen in Fig. Ul we 
find a steady decrease in the value of d as iV increases, 
that is, our distribution approaches a Gaussian. 

Finally, we turn to the velocity amplitude, which while 
extremely well reproduced in the effective undccoratcd 
chain is not captured well by the binary collision approx- 
imation. It is easy to see the source of the problem. 
Recall that Eq. (fT9)) is valid only when the masses of the 
granules in the effective chain are assumed to all be equal 
and only the interactions between them are affected by 
the randomness of the sizes of the smaller granules. The 
agreement between the theory, Eq. (|19[) . and the numer- 
ics indicates that these assumptions are valid when calcu- 
lating pulse travel times. However, the velocity obtained 
from the model as posed in Eq. (|14[) is independent of k 
(and only very weakly dependent on k if the direct < 6% 
effect of the randomness of the masses is not neglected) 
and does not follow the behavior observed in the numer- 
ical solution of the exact equations, Fig. [2] The problem 
lies in the fact that the velocity Vk obtained from the 
binary model depends only on the mass ratio of two col- 
liding granules. However, the random interactions can 
introduce unavoidable three (or more) granule scenarios. 
For example, if the interaction between a granule 1 and a 
granule 2 is weak but that between 2 and 3 is strong, then 
a collision between 1 and 2 will inevitably involve gran- 
ule 3. A remedy might be to include this effect through 
a further rcnormalization of the masses resulting from 
the random interactions. This will be explored in future 
work. 



V. CONCLUSIONS 

We have studied pulse propagation in ID granular 
chains decorated with small granules of random radii in- 
serted between each pair of large granules all of the same 
size. This study has proceded in two steps. 

Firstly, we used the effective scheme introduced in 
Ref. [31( to obtain an equivalent undecorated random 
chain, and showed that this effective description works 
remarkably well for all properties tested, so well that the 
numerical results obtained from the original chain and 
from the effective chain are essentially indistinguishable. 
Since the effective chain is only half as long as the original 
chain, this represents a considerable savings in computa- 
tional effort. 

Secondly, using the binary collision approximation on 
the effective chain, we have obtained analytic expressions 
for the velocity amplitude of the pulse and the time that 
the pulse takes to reach the fcth granule (pulse speed 
dk/dt) along the chain. By construction, the analytic re- 
sults obtained using the binary collision approximation 
neglect the effects of the randomness of the small gran- 
ules on the masses of the large granules in the effective 
chain. The randomness appears only in the interactions 
between granules. It is thus not surprising that the ve- 
locity amplitude of the pulse obtained from the binary 
model does not predict the correct behavior as seen in nu- 
merical results, since this amplitude is determined by the 
masses of colliding granules and is therefore affected by 
the randomness of these masses. In the previous section 
we have suggested a possible remedy to this issue. On 
the other hand, the pulse speed and the distribution of 
the times taken by the pulse to reach the fcth granule are 
very well reproduced by the theory. This is because the 
time of pulse propagation depends on the interactions, 
being shorter (longer) for stronger (weaker) interactions. 
Our theory incorporates this dependence very accurately, 
cf. Eq. JT5J. 

As was noted in Ref. [3l|, the effective description 
works well as long as the size of the small granule re- 
mains less than ~ 40% of the bigger granule. This places 
a restriction on the size/mass distribution of the smaller 
granules. Here, for simplicity, we have considered a uni- 
form distribution of the smaller granules. However, the 
validity of the effective description and of the binary- 
collision approximation will remain valid for arbitrary 
distributions as long as the size restriction on the small 
granules is satisfied. In future work we plan to generalize 
our effective chain description (and the associated binary 
collision approximation) to other chain configurations, 
with the eventual goal of understanding pulse propaga- 
tion in chains of arbitrary granular configurations. 
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